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1.0  INTRODUCTION 

The  interference  flow  field  about  a  store-laden  aircraft  is  frequently  the  most  important 
parameter  affecting  a  store  trajectory.  Lack  of  accurate  experimental  data  has  been  a  major 
hindrance  to  development  and  evaluation  of  analytical  techniques  for  computing  such 
flows,  especially  in  the  transonic  regime.  The  occurence  of  shock  waves  within  the  flow 
field,  often  accompanied  by  regions  of  flow  instability,  makes  measurement  by  conventional 
mechanical  probes  unreliable.  A  nonperturbing  device  such  as  a  laser  Doppler  velocimeter 
(LV)  must  be  used  if  accurate  velocity  measurements  are  to  be  obtained  throughout  a 
transonic  flow  field. 

The  purpose  of  this  investigation  was  to  obtain  LV  measurements  in  the  flow  fields  of 
several  wing/ pylon/store  configurations.  The  flow-field  surveys  were  conducted  with  use  of 
a  three-component  LV  system  in  the  Arnold  Engineering  Development  Center  (AEDC) 
Aerodynamic  Wind  Tunnel  (IT)  of  the  Propulsion  Wind  Tunnel  (PWT)  Facility. 

2.0  APPARATUS 


2.1  TEST  FACILITY 

The  AEDC  Tunnel  IT  is  a  continuous-flow,  open-circuit  wind  tunnel  that  can  be 
operated  over  a  Mach  number  range  from  0.20  to  1.50.  The  tunnel  operates  at  a  constant 
stagnation  pressure  of  approximately  2,850  psfa.  Stagnation  temperature  in  the  tunnel  can 
be  varied  from  80  to  120°F  above  ambient  temperature.  A  complete  description  of  the 
facility  is  included  in  Ref.  1. 

The  four  test  section  walls  used  during  the  test  were  6-percent  perforated  except  for  a 
closed  region  on  either  sidewall  extending  from  Sta  14.6  to  Sta  26.6.  Optical-quality 
windows  were  installed  in  the  plenum  and  test  section  sidewalls  adjacent  to  the  LV.  A 
variable  angle-of-attack  wing  mount  was  located  between  the  same  tunnel  stations  on  the 
opposite  test  section  wall.  The  plenum  window  was  mounted  in  a  special  inset  assembly  that 
allowed  positioning  it  close  to  the  tunnel  window.  The  arrangement  permitted  the  LV  beam 
array  to  be  traversed  about  the  test  article  without  requiring  a  large,  prohibitively  expensive 
plenum  window.  A  plan  view  of  the  Tunnel  IT  test  area  showing  the  LV  system  installation 
is  presented  in  Fig.  1 ,  The  fluidized  bed  particle  seeder  shown  exhausting  into  the  air  intake 
room  was  utilized  to  supplement  the  ambient  atmospheric  particles  and,  thus,  to  improve 
the  data  acquisition  rate. 
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2.2  TEST  ARTICLES 

Test  hardware  included  a  wall-mounted  45-deg  swept  wing  with  a  constant  6-in.  chord 
NACA  0005-34  airfoil  cross  section,  two  1/20-scale  store  models  (an  MK-83  and  an  M-l  17), 
and  two  pylons,  one  with  unswept  and  the  other  with  60-deg  swept  leading  and  trailing 
edges.  Model  dimensions  are  presented  in  Fig.  2.  The  configuration  identification  and 
coordinate  system  are  given  in  Figs.  3  and  4,  respectively. 

Fourteen  test  configurations  were  investigated.  The  configurations,  defined  in  Table  1, 
are  all  variations  of  the  five  basic  wing/pylon/store  combinations  shown  in  Fig.  5. 

2.3  INSTRUMENTATION 

The  standard  tunnel  system  was  used  to  record  the  test  conditions  and  model  pressures. 
The  flow-field  measurement  utilized  a  three-velocity  component  LV  system.  Figure  6  shows 
the  optical  arrangement  of  the  LV.  The  system  consisted  basically  of  both  a  two-component 
LV  and  a  single-component  LV  mounted  on  the  same  traverse  table  and  sharing  a  common 
laser  light  source.  A  blue  beam,  two-component  LV  provided  the  magnitude  of  the  velocity 
components  parallel  to  the  x  and  z  axes  of  the  mutually  orthogonal  wind  tunnel  coordinate 
system.  Rather  than  measuring  the  component  magnitude  parallel  to  the  y  axis  directly,  a 
green  beam,  single-component  LV  obtained  the  velocity  component  in  the  y'  direction.  A 
more  detailed  description  of  the  system  is  contained  in  Appendix  A.  The  technique  for 
resolving  the  y  component  of  velocity  from  the  measured  components  is  described  in 
Appendix  B  and  Ref.  3. 

2.4  FLUIDIZED  BED  PARTICLE  SEEDER 

The  tunnel  flow  was  seeded  with  alumina  particles  with  a  nominal  diameter  of  one 
micron.  As  shown  in  Fig.  1 ,  the  seed  particles  were  introduced  into  the  air  intake  room  along 
with  atmospheric  particles.  These  particles  then  passed  through  the  intake  filters  and  the 
compressor,  which  tended  to  remove  the  larger  particles  and  break  up  clustered  seed 
particles.  Under  unseeded  conditions  the  time  required  to  record  an  LV  data  point  varied 
from  several  seconds  to  several  tens  of  seconds  depending  upon  the  atmospheric  dust 
content.  Operation  of  the  seeder  reduced  this  time  by  approximately  50  percent.  It  is 
important  to  emphasize  that  the  amount  of  seed  material  required  was  minimal:  Less  than 
three  pounds  of  alumina  powder  were  used  during  several  hundred  hours  of  testing.  Neither 
the  tunnel  nor  model,  apparently,  suffered  any  adverse  effects. 
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3.0  PROCEDURE 


3.1  TEST  CONDITIONS 

The  test  was  designed  to  investigate  the  mutual  interference  between  components  of  each 
configuration  in  the  transonic  Mach  number  range.  As  shown  in  Table  1  each  configuration 
was  tested  at  M»  =  0.920,  which  was  determined  experimentally  (Ref.  2)  as  the  lowest  Mach 
number  for  which  a  shock  was  firmly  established  between  the  store  and  the  wing  lower 
surface  of  Configuration  201.  In  addition,  selected  configurations  were  tested  at  M„  = 
0.975  and/or  1.025.  During  the  tests,  a  tunnel  stagnation  temperature  of  160°F  was 
maintained  to  minimize  velocity  continuity  errors. 

3.2  DATA  ACQUISITION 

The  LV  data  were  recorded  by  a  computer-based  data  acquisition  system  shown  in  block 
diagram  form  in  Fig.  7.  Nominally  1,000  velocity  samples  were  recorded  for  each  of  the 
velocity  components  at  each  point  of  interest  in  the  flow  field.  The  time  required  to  record 
these  data  varied  from  several  seconds  to  several  tens  of  seconds,  depending  upon  the 
atmospheric  dust  content  of  the  intake  air.  A  more  desirable  and  controllable  rate  of  data 
acquisition  was  obtained  by  using  the  seeder. 

The  terminal  shown  in  Fig.  7  provided  an  online  view  of  the  LV  velocity  data  and  the 
free-stream  tunnel  velocity  derived  from  tunnel  pressure  and  temperature  measurements. 
The  printer  provided  a  hard  copy  of  these  CRT  data  while  the  magnetic  tape  unit  was  the 
primary  means  for  recording  data. 

The  focal  volume  position  relative  to  the  test  model  was  recorded  for  each  velocity  data 
set.  Periodic  corrections  were  necessary,  however,  to  minimize  spatial  error  caused  by 
temperature  effects. 

3.3  DATA  REDUCTION 

A  generalized  LV  data  reduction  and  analysis  software  package  has  been  developed  that 
can  process  up  to  three  components  of  either  nonsimultaneous  or  simultaneous  velocity 
data.  The  program  provides  for  extensive  statistical  analysis,  including  computations  of 
sample  statistics  and  estimations  of  population  parameter  values.  The  quantities  computed 
include  mean,  variance,  standard  deviation,  third  and  fourth  moment,  skewness,  kurtosis, 
and  absolute  deviation  as  defined  in  Ref.  4.  The  mode  of  each  sample  distribution  can  also 
be  determined  independent  of  the  performance  of  “split”  statistical  computations  on  the 
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portions  of  the  distribution  on  either  side  of  the  mode  value.  For  nonsimultaneous  data,  the 
component  means  and  modes  can  be  transformed  from  LV  measurement  coordinates  into 
any  desired  mutually  orthogonal  axis  system  and  the  total  mean  and  mode  velocity  vector 
magnitudes  can  also  be  calculated  (Appendix  B).  For  simultaneous  data  sets,  the  coordinate 
transformation  can  be  performed  on  the  discrete  particle  velocity  component  sets  before  the 
statistical  analyses.  The  correlation  provided  by  the  simultaneity  condition  also  permits  the 
computation  of  several  additional  statistical  parameters  including  covariances  (Reynolds 
shear  stresses),  correlation  coefficients,  the  generalized  variance,  and  the  scatter  coefficient 
(Ref.  4). 

A  master  summary  file  containing  the  results  of  the  various  data  reduction  procedures 
can  be  accessed  for  priming  or  plotting  output  or  for  further  analysis.  A  diagram  of  the 
offline  data  processing  sequence  is  presented  in  Fig.  8. 

3.4  PRECISION  OF  MEASUREMENTS 

3.4.1  Tunnel  Conditions 

Care  was  taken  to  allow  the  tunnel  to  stabilize  before  the  relatively  long  data-recording 
period.  The  resulting  variation  in  mean  free-stream  Mach  number  was  estimated  to  be 
±  0.003.  The  variation  in  mean  free-stream  velocity  was  less  than  ±  5.0  ft/sec.  The  sidewalls 
with  windows  presented  an  additional  source  of  Mach  number  uncertainty.  According  to  the 
Tunnel  IT  calibration  (Ref.  5),  use  of  the  windows  can  result  in  significant  variation  in  Mach 
number  along  the  tunnel  centerline.  A  statistical  analysis  of  the  centerline  Mach  number 
resulted  in  a  2 -a  deviation  from  the  mean  free-stream  Mach  number  of  ±0.035. 

3.4.2  Spatial  Resolution 

The  LV  probe  volume  was  positioned  in  the  test  article  flow  field  by  moving,  as  a  unit, 
the  laser  and  optical  components  shown  on  the  traverse  table  in  Fig.  6.  The  traverse  system 
could  be  programmed  to  step  in  integer  multiples  of  0.001  in.  in  the  x-  and  y-axis  directions 
and  in  integer  multiples  of  0.00038  in.  in  the  z-axis  direction.  The  tip  of  the  store  model 
served  as  the  origin  of  the  x-y-z  tunnel  reference  coordinates.  Because  of  thermal  expansion 
effects,  this  origin  moved  as  a  function  of  elapsed  tunnel  run  time.  Therefore,  it  was 
necessary  to  relocate  the  model  tip  periodically  by  observing  the  precisely  repeatable 
diffraction  patterns  by  the  transmitted  laser  beams,  contacting  the  tip  and  upper  and  lower 
surfaces  of  the  store  model.  The  repeatability  of  the  model  tip  position  determinations  was 
within  ±0.010  in.  for  the  x  and  z  axes. 


8 


A  EDC-TR-80-54 


A  velocity  scan  was  the  means  to  determine  the  lateral  (y)  position  of  the  model.  Laser 
velocity  data  were  recorded  at  from  IS  to  20  points  along  a  constant  y-z  line  from  one  side  of 
the  model  to  the  other  in  the  region  of  high-velocity  gradient  near  the  model  nose.  The 
actual  model  centerline  position  could  usually  be  determined  to  within  ±  0.030  in.  from  the 
symmetry  of  the  velocity  distribution  in  the  y  direction.  The  technique  is  suitable  for  most 
models  symmetrical  about  a  constant  y  plane.  The  results  are  valid  at  any  angle  of  attack  as 
long  as  the  yaw  angle  remains  small;  this  angle  is  critical  because  an  error  in  yaw  caused  by 
tunnel  flow  angularity  or  model  misalignment  can  destroy  the  symmetry  of  the  upward 
velocity  distribution,  resulting  in  a  larger  y-position  uncertainty.  It  should  be  noted  that  all 
the  Vy  distributions  observed  in  the  present  test  were  acceptably  symmetrical. 

3.4.3  Laser  Velocimeter  Measurements 

The  velocity  measurements  in  thex,  y '  ,  .and  z  directions  were  usually  repeatable  to  within 
±10.0  ft/sec.  This  corresponds  to  the  relative  uncertainties  of  ±1.0  and  ±  1 .2  percent  in  the 
Vux  and  Vy  components,  respectively,  at  M«,  =  0.92.  The  uncertainty  in  the  resolved  Vy 
component  was  found  to  be  somewhat  higher  than  that  observed  in  the  components 
measured  directly.  The  error  sources  are  discussed  in  Appendix  C. 

4.0  RESULTS  AND  DISCUSSION 

LV  velocity  measurements  were  made  for  each  test  configuration  at  the  same 
predetermined  survey  points  in  the  vicinity  of  the  test  article.  The  particular  points  were 
selected  to  provide  optimum  flow-field  definition  and  to  assess  mutual  aerodynamic 
interference  between  the  model  elements.  Several  of  the  survey  grid  points  were  located  on 
the  opposite  side  of  the  pylon  from  the  LV  system  and  required  measurement  either  through 
the  transparent  pylon  or  through  the  gap  between  the  pylon  and  store  (Fig.  9).  All  survey 
point  coordinates  are  referenced  to  the  store  nose  as  illustrated  in  Fig.  4. 

Velocity  vector  projections  measured  in  several  constant  x  planes  are  presented  in  Fig.  10 
for  Configuration  111.  Figures  10c  and  d  contain  examples  of  measurements  taken  through 
the  pylon  and  measurements  obtained  through  the  pylon/store  gap. 

The  vector  projections  of  velocities  measured  in  a  constant  z  plane  passing  through  the 
pylon/store  gap  are  presented  in  Figs.  11  through  13.  Figures  11  and  12  show  that  the 
unswept  pylon  causes  more  severe  flow-field  disturbances  in  the  forward  region  than  does 
the  swept  pylon;  the  figures  also  show  that  larger  sidewash  angles  occur  beneath  the  aft 
portion  of  the  pylon  for  a  wing  angle  of  attack  of  zero  (Configurations  121  and  111).  The 
sidewash  disturbance  is  stronger  at  a  wing  angle  of  attack  of  2.0  deg,  however,  for  the 
unswept  pylon  (Configuration  221). 
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Mach  number  effects  on  the  Configuration  121  flow  field  can  be  seen  in  Figs.  12b  and  13. 
Although  stronger  sidewash  appears  near  the  pylon  at  M*.  =  0.92  than  at  M#  =  0.975,  the 
subsonic  flow  fields  are  quite  similar  at  z/D  =  0.56.  The  supersonic  case,  shown  in  Fig.  13b, 
produced  a  general  outward  flow  tending  toward  the  wing  tip. 

The  effect  of  pylons  upon  flow-field  velocity  vectors  measured  in  a  constant  y  plane  are 
presented  in  Figs.  14  and  15  for  the  MK-83  and  M-l  17  store,  respectively,  with  the  wing  at 
zero  deg  angle  of  attack.  Except  for  the  differences  in  vector  angle  and  magnitude  near  the 
pylon  leading  edge,  the  MK-83  flow  fields  appear  to  be  nearly  identical.  Likewise,  the  M-l  17 
configurations  exhibit  only  minimal  differences  in  the  selected  constant  y  plane.  The 
gradients  are  generally  larger  for  the  M-l  17  than  for  the  MK-83  configurations. 

Distributions  of  the  three  orthogonal  components  of  velocity  measured  “along  a  line 
where  the  centerline  of  the  store  would  be  if  it  were  installed”  (Ref.  6)  are  presented  in  Fig. 
16  for  the  three  configurations  with  a  wing  angle  of  attack  of  2.0  deg.  Such  data  could  be 
used  for  computing  initial  or  launch  conditions  for  a  store  using  the  flow  angularity 
techniques  described  in  Refs.  6  and  7.  As  previously  observed,  the  unswept  pylon  generates 
larger  flow  field  disturbances  than  does  the  swept  pylon. 

In  addition  to  obtaining  the  flow-field  grid  data,  more  detailed  surveys  were  made  for 
selected  configurations  to  detect  the  presence  and  location  of  embedded  shock  waves.  The 
results  of  three  such  surveys  (Fig.  17)  illustrate  the  effect  of  wing  angle  of  attack  on  the 
shock  characteristics  of  Configurations  101  and  201.  The  velocity  distributions  at  z/D  = 
0.74  reveal  what  appears  to  be  a  secondary  weak  shock  at  approximately  z/D  =  4.0  in  the 
Configuration  101  flow  field;  this  shock  is  not  present  for  Configuration  201.  A  survey 
closer  to  the  store  at  z/D  =  0.56  confirmed  the  existence  of  one  and  possibly  two  additional 
shock-like  disturbances  forward  of  the  strong  primary  shock  at  x/D  =  5.7.  The  velocity 
distribution  immediately  downstream  of  the  Configuration  201  shock  is  essentially  the  same 
as  that  observed  previously  in  the  two-component  LV  investigation  (Ref.  2),  indicating  that 
similar  errors  from  particle  dynamics  were  present.  Discussions  of  the  so-called  “particle 
lag”  effect  can  be  found  in  Refs.  2  and  8. 

The  LV  boundary-layer  profiles  described  in  Fig.  18  were  obtained  along  a  constant  x-y 
line  beneath  the  wing  of  Configuration  100  at  the  40-percent  chord  point.  The  results  are 
plotted  against  the  local  normal  coordinate,  zs/S,  which  is  positive  away  from  the  surface 
and  nondimensionalized  by  an  estimated  boundary-layer  thickness  of  0.16  in.  The  large 
values  of  (y/VE)  and  [(VxVy)ME/VE2]  which  occur  in  the  profiles  at  (Zj/6)  >  1.0  are  not 
consistent  with  the  assumption  of  isotropic  turbulence  outside  the  boundary  layer;  these 
large  values  are  probably  caused  by  the  y-component  uncertainty  discussed  in  Appendix  C. 
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If  the  higher-order  statistical  information  presented  in  Figs.  18b  and  18c  is  of  primary 
interest,  the  LV  measurement  coordinate  system  should  be  optimized  to  minimize 
y-component  resolution  uncertainty. 

5.0  CONCLUDING  REMARKS 

The  transonic  flow  fields  of  several  wing/pylon/store  configurations  were  investigated 
by  using  a  three-velocity-component  LV.  The  results  indicate  that  an  LV  can  provide  useful 
nonintrusive  velocity  measurements  for  such  applications. 

The  representative  examples  of  the  experimental  data  presented  reveal  the  flow  field  to 
be  rather  complex  in  constrast  to  the  relative  simplicity  of  the  model.  As  a  rule,  the  unswept 
pylon  generated  larger  disturbances  than  did  the  swept  one.  Larger  flow  deflections  were 
observed  for  the  M-l  17  store  than  were  observed  for  the  less  blunt  MK-83.  The  disturbance 
that  occurred  near  the  lower  pylon  leading  edge  was  intensified  by  the  presence  of  the  store. 

The  present  investigation  marked  the  first  use  of  a  three-component  LV  at  the  AEDC. 
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Figure  1.  Tunnel  IT  layout  and  LV  system  installation. 
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Figure  2.  Model  dimensions. 
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Figure  2.  Concluded. 
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CONFIGURATION 


W  =  0  NO  WING 

1  WING  AT  aw=  0  deg 

2  WING  AT  aw  =  2.0  deg 

P  =  0  NO  PYLON 

1  UNSWEPT  PYLON 

2  SWEPT  PYLON 

S  =*  0  NO  STORE 

1  MK  -  83 

2  M  -  117 

Configuration  convention  of  Ref.  2  was  modified 
to  include  pylons  and  incidence  angles. 

Figure  3.  Configuration  identification. 
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PLAN  VIEW 
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a.  Configuration  221 


b.  Configuration  211 


d.  Configuration  111 


e.  Configuration  122 

Figure  5.  Basic  wing/pylon/store  configurations. 
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Figure  6.  Tunnel  IT  three-velocity  component  LV  optical  system  installation. 
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Figure  7.  LV  data  acquisition  system  schematic. 
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Figure  8.  Offline  data  processing  sequence. 
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Figure  9.  Configuration  221  installed  in  tunnel  IT  with  the  three-component  LV  beams 
passing  through  the  gap  between  the  pylon  and  store. 
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The  crosshatched  area  Indicates  model  cross  section  at  this  x/D  station. 


a.  x/D  =  0.20 

Figure  10.  Mean  velocity  vector  projections  in  constant  x  planes  for 
configuration  111  <aw  =  0),  =  0.92. 
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a.  Configuration  111 

Figure  12.  Effects  of  pylon  on  velocity  vector  projections  in  constant  z  plane  for 
four  wing/MK-83  configurations  {a*,  =0),  =  0.92,  z/D  =  0.56. 
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b.  Configuration  121 
Figure  12.  Continued. 
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Figure  13.  Effect  of  Mach  number  on  mean  velocity  vector  projections  in  constant  z  plane 
for  three  wing/MK-83  configurations  (aw  =  0),  =  0.92,  y/D  =  0. 


SCALE 


1 - 1 - 1 _ I _ I 

0  2  4  6  8 

x/D 


a.  Configuration  111 

Figure  14.  Effect  of  pylon  on  velocity  vector  projections  in  constant  y  plane  for 
three  wing/M  K-83  configurations  {a*  =  0),  =  0.92,  y/D  =  0. 
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b.  Configuration  121 
Figure  14.  Continued. 


AEDC-TR-aO-64 


z/D 


SCALE 
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Figure  15.  Effect  of  pylon  on  velocity  vector  projections  in  constant  y  plane  for 
two  wing/M-117  configurations  (a*  =  2.0),  =  0.92,  y/D  =  0. 
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Figure  17.  Effect  of  wing  angle  on  the  distribution  of  mode  velocity  magnitude 
between  the  wing  and  store,  =  0.92,  y/D  =  0. 
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b.  Turbulent  intensity 

Figure  18.  Wing  boundary-layer  characteristics  on  Configuration  100  lower  surface, 
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Appendix  A 

Laser  Doppler  Velocimeter  System  Details 

The  LV  shown  in  Fig.  6  is  a  differential  Doppler  system  which  is  further  classified  as  a 
moving  fringe  system  operating  in  the  backscatter  mode.  The  laser  was  operated  in  the  all¬ 
lines  mode  at  from  6  to  8  w,  although  the  maximum  power  capability  was  18  w.  In  the  all¬ 
lines  mode  the  output  laser  beam  actually  consists  of  nine  different  wavelengths. 

The  blue  (0.488-jun  wavelength)  and  green  (0.5145-/tm  wavelength)  beams  are  the  two 
most  powerful  of  the  nine.  These  two  beams  were  separated  from  the  rest  by  color  filters 
CF1  and  CF2  (Fig.  6).  The  energy  of  the  seven  remaining  beams  was  dissipated  by  the  beam 
dump,  BD.  The  blue  beam  was  passed  through  an  acousto-optical  device  termed  a  dual 
Bragg  cell  (Ref.  9),  represented  in  Fig.  6  as  the  element  DBC.  The  DBC  splits  the  single  input 
beam  into  four  nominally  equal  power  output  beams  as  shown  by  beam  pattern  3,  4,  5,  6  in 
Fig.  6.  Nominally  25  percent  of  the  input  beam  passes  through  the  DBC  and  retains  its 
original  optical  frequency,  v,  where  v  =  C/X,  C  is  the  speed  of  light,  and  X  is  the  wavelength 
of  the  laser  light.  Interaction  of  the  input  laser  beam  with  acoustic  waves  in  the  DBC  causes 
beams  4,  5,  and  6  to  be  split  off  from  beam  3.  Interaction  with  the  acoustic  waves  also  causes 
the  optical  frequency  of  beam  4  to  shift  to  v  +  15  MHz,  beam  5  to  v  +  60  MHz,  and  beam  6 
to  v  +  45  MHz.  These  four  beams  are  projected  into  the  tunnel  test  section  and  caused  to 
intersect  by  lens  L9.  Lenses  LI,  L2,  and  L3  cause  the  green  and  blue  beams  to  focus  in  the 
intersection  region.  The  diameter  of  their  cigar-shaped  intersection  regions  was  about  300 
jim  for  both  the  blue  and  green  beams. 

Beam  pairs  3,4  and  5,6  interfere  in  the  intersection  region  and  produce  alternate  bright 
and  dark  fringe  planes  that  move  in  the  x  direction  at  a  15-MHz  rate,  which  is  the  frequency 
of  the  electrical  drive  signal  to  the  x-axis  Bragg  cell  transducer.  A  stationary  particle  in  the 
intersection  region  scatters  light  at  a  15-MHz  rate  as  the  bright  fringes  sweep  past  it.  A 
particle  traveling  in  the  +x  direction  scatters  light  at  a  rate  greater  than  15  MHz,  whereas  a 
particle  traveling  in  the  -x  direction  scatters  light  at  a  rate  less  than  15  MHz.  The  moving 
fringes,  therefore,  resolve  the  directional  ambiguity  problem  associated  with  static  fringe  LV 
systems. 

The  z-axis  measurement  capability  is  furnished  by  beam  pairs  3,6  and  4,5,  which  produce 
fringes  that  travel  in  the  -z  direction  at  a  rate  of  45  MHz.  Beam  pairs  3,5  and  4,6  also 
interfere  and  thereby  define  diagonal  axes  with  fringes  moving  at  rates  of  30  and  60  MHz, 
but  the  signals  derived  from  these  two  axes  are  later  rejected  by  electronic  filtering 
techniques.  Scattered  light  is  collected  by  lenses  L5  and  L6  and  focused  at  aperture  API 
Color  filter  CF3  passes  the  blue  and  rejects  the  green  laser  light.  The  photomultiplier  tube 


47 


A  E  DC-TR-80-54 


PMT1  then  produces  an  electrical  signal  in  response  to  the  scattered  blue  laser  light.  The 
x-axis  signal  centered  at  15  MHz  and  the  z-axis  signal  centered  at  45  MHz  are  separated  by 
bandpass  electronic  filters  (Ref.  10).  Electronic  signal  processing  equipment  (Ref.  1 1)  is  then 
used  to  determine  the  frequencies  of  the  x-axis  and  z-axis  signals  that  are  proportional  to  the 
particle  vector  component  velocities. 

The  description  of  the  green  beam  array  is  similar  to  that  of  the  blue  except  that  only  one 
Bragg  cell  is  needed  to  produce  the  two  green  beams.  Motion  of  the  green  fringes  is  at  a 
15-MHz  rate  in  the  -y'  direction.  Scattered  light  is  collected  by  lenses  L7  and  L8,  and  the 
green  light  is  detected  by  PMT2.  The  x-y-z  velocity  components  are  then  derived  by 
transformation  from  the  x-y'-z  axis  measurements. 

The  angle  subtended  by  the  axes  of  the  transmitted  blue  and  green  beams  was  set  at 
nominally  30  deg,  and  the  angle  subtended  by  a  receiver  axis  and  its  associated  transmitter 
axis  was  set  at  nominally  20  deg.  The  intersection  of  the  transmitter  focal  volumes  and  the 
receiver  image  volumes  defines  the  spatial  volume  from  which  the  scattered  light  signals 
originate.  This  volume  is  called  the  probe  volume.  The  dimension  of  the  probe  volume  in  the 
x  and  z  directions  was  approximately  300  pm  and  the  dimension  in  the  y  direction  was 
approximately  600  pm. 
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Appendix  B 

Resolution  of  the  y  Component  of  Velocity 

The  signal  frequency  produced  by  one  axis  of  a  differential  Doppler-type  LV  system  is 
given  by 


f 


D 


(B-l) 


where  V  is  the  velocity  vector  of  the  light-scattering  particle,  and  k|  and  kj  are  the 
propagation  vectors  of  the  two  intersecting  laser  beams.  From  Eq.  (B-l)  it  can  be  readily 
seen  that  an  LV  system  simply  measures  the  orthogonal  projection  of  V  onto  an  axis  defined 
by  k2  -  k|.  As  with  the  present  case,  it  is  oftentimes  technically  impractical  to  implement  a 
three-axis  LV  system  in  which  the  three  axes  are  mutually  orthogonal.  In  general,  the 
relationship  of  the  vector  components,  V;,  of  V  referred  to  the  orthogonal,  Cartesian,  tunnel 
reference  coordinates,  Xj,  in  terms  of  the  (V  ±)j-  projections  of  V  onto  the  oblique,  Cartesian 
coordinates,  V ,  defined  by  the  LV  measurement  axes  i' ,  given  by  (Ref.  3) 


V 

(C2'2C3'3  c3  '2C2  '3)  (C3'2C1'3  c]'2c3'a)  (cr2c2 '3  c2  *2C  1  '3  ]) 

'(v4>; 

V2 

1 

“  7\ 

(c3'lc2'3  “  c2'lc3  '3)  (C1  '  1  c3  '3  ~  c3  '1 C1  '3)  (c2'lcl'3  “  C1  'l  c2  '3  ) 

('1)2' 

J1  3_ 

(^•lc3'2  ~  c3'lc2  '2)  (C3'1C1'2  "  C1  '1  c3  '2)  (CI'IC2'2  ~  C2'lcr2)_ 

.(VJ  3 : 

(B-2) 

where 


-N  =  cr]  (C2'2C3  '1  W3'2)  '  ^’2  (^'^'3 


3)  +  c[\3  (c2'l 


3  12  2  J 


(B-3) 


and  Cj'j  is  the  cosine  of  the  angle  defined  by  the  intersection  of  the  x;'  and  Xj  axes.  In  the 
present  case 


c 


l'l 


~3  '3 


=  1 


'(B-4) 


and 

C1'2  =  C1'3  =  c2'3  =  c3'l  =  c3'2  =  0  (B-5) 

since  the  X| '  and  X3'  LV  axes  are  coincident  with  the  X|  and  X3  tunnel  reference  axes.  Now 
with  more  familiar  notation  where  Xj  =  x,  x2  =  y,  x3  =  z,  and  xj '  =  x' ,  x2'  =  y',  X3'  = 
z',  Eq.  (B-2)  becomes 
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Appendix  C 

Uncertainty  in  y  Component  of  Velocity 


AEDC-TR-80-64 


The  overall  uncertainty,  Uf,  in  the  determination  of  a  quantity  f(pj)  that  is  functionally 
dependent  upon  various  parameters,  pj,  can  be  expressed  by  (Refs.  12  and  13) 

Of  -  1  (Dl  *  ‘0.95  Sf)  (C-l) 


where  Bf  is  the  overall  bias  error,  Sf  is  the  overall  precision  error,  and  to.95  is  the  95th 
percentile  point  for  the  Student  t-distribution.  For  the  degrees  of  freedom  typical  of  an  LV 
data  set,  to.95  *  2.0,  the  overall  bias  error  is  defined  by 


(C-2) 


where  bpj  is  an  elemental  bias  error.  Similarly,  the  overall  prevision  error  is  related  to  the 
elemental  precision  errors,  sK,  by 


(C-3) 


The  overall  uncertainty  and  repeatability  in  measurements  of  a  single  velocity  component 
obtained  with  the  LV  equipment  used  at  AEDC  have  been  previously  determined  (Refs.  2 
and  14).  The  measurement  uncertainty  and  repeatability  in  measurement  of  a  single  particle 
velocity  in  the  transonic  velocity  range  should  be  within  ±  1.5  percent  of  the  recorded  value 
while  the  uncertainty  and  repeatability  of  mean  velocity  measurements  should  be  within 
±1.3  percent.  However,  measurement  repeatability  fluctuations  of  several  percent  were 
noted  whenever  the  three-component  LV  system  was  subjected  to  varying  ambient 
temperature  conditions.  Temperature-change-induced  movements  of  the  optical  elements 
preceding  the  two  Bragg  cells  caused  changes  in  the  spatial  orientation  of  the  beams  entering 
the  Bragg  cells,  ultimately  causing  small  changes  in  the  probe  volume  fringe  spacings 
(Appendix  A).  The  intensity  of  the  incident  laser  beams  and  the  low  mass  and  low  thermal 
conductivity  of  the  color  filters,  CF1  and  CF2  (Fig.  6),  probably  accounted  for  most  of  the 
observed  temperature  effects.  The  measurement  fluctuations  typically  subsided  to  the 
±  1 -percent  range  after  the  tunnel  and  surrounding  area  reached  a  stable  operating 
temperature  as  long  as  gross  air  currents  in  the  vicinity  of  the  LV  system  were  suppressed. 


With  use  of  Eqs.  (C-l),  (C-2),  and  (C-3),  an  expression  for  the  uncertainty  in  the 
determination  of  the  transformed  component  Vy  can  be  developed  as  (with  the  assumption 
of  uncorrelated  elemental  errors) 
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(‘■7;)‘W.]' 4  [teh'.).-] 

-v.l 


(C-4) 


As  might  be  expected,  this  expression  shows  that  the  uncertainty  is  a  minimum  when  the  y' 
axis  is  coincident  with  the  y  axis  and  is  a  maximum  when  the  y'  axis  is  coincident  with  the  x 
axis.  For  the  three-component  LV  system,  the  nominal  angle  subtended  by  y '  and  y  was  120 
deg,  whereas  the  angle  subtended  by  y '  and  x  was  30  deg  with  an  uncertainty  in  both  angles 
of  ±0.3  deg.  With  use  of  previously  stated  uncertainties  for  (Vx)x'  and  (Vj_)y\  Eq.  (C-4) 
for  the  uncertainty  in  the  determination  of  Vy  for  a  single  particle  from  simultaneous  (V ± )x ' 
and  (Vx)y-  measurements  is 

«».  -  ij[(o.026)(vi)I.]2  +  [(-0.003)  (vj,-]2  4  [(o.00006)(v±)„.]2 

)  (C-S) 

:  [(0.010) (vj,.-  (O.OIZ) 

For  small  values  of  (V  x)y ' ,  the  first  and  last  terms  are  dominant  and  impart  an  uncertainty 
of  about  +(0.036)  (VXV  to  determinations  of  Vy.  Thus,  the  uncertainty  amplifying  the 
factor  (Cy  ’  x)/(Cy '  y)  becomes  a  primary  concern  when  small  values  of  (Vx)y'  are  to  be 
measured  and  when  the  angles  subtended  by  the  y'  and  y  axes  diverge  from  0  or  180  deg  and 
approach  90  or  270  deg,  respectively. 

Fortunately,  because  the  elemental  errors  are  uncorrelated,  the  error  in  mean  velocity 
decreases  with  increasing  sample  size.  The  uncertainty  does,  however,  affect  the  higher- 
order  statistical  parameters  since  the  sample  histogram  is  artificially  broadened. 
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NOMENCLATURE 

Axial  distance  from  store  nose,  in.  (Fig.  2a) 

Overall  bias  error  in  determining  f(pj) 

Elemental  bias  error 
Speed  of  light 

Cosine  of  the  angle  between  the  Xj '  and  Xj  axes 
Maximum  store  diameter,  in. 

Local  store  diameter,  in.  (Fig.  2a) 

Pylon  contour  coordinate,  in.  (Fig.  2d) 

Pylon  contour  coordinate,  in.  (Fig.  2d) 

Doppler  frequency 

Function  dependent  upon  parameters,  p; 

Propagation  vector  of  the  ith  laser  beam 

Free-stream  Mach  number 

Overall  precision  error  in  determining  f(pj) 

Elemental  precision  error 
Tunnel  station,  in. 

Half-thickness  of  wing,  in.  (Fig.  2c) 

Uncertainty  in  the  determination  of  the  function,  f(pj) 
Velocity  vector 

Velocity  vector  magnitude,  ft/sec. 

Velocity  at  boundary-layer  edge,  ft/sec. 

; _ 

ith  component  of  V,  ft/sec. 

Projection  of  V  onto  Xi'  axis,  ft/sec 
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(Vi)x' 

(ViV 

(VOz' 

Vx,  Vy,  vz 
v;,  v;,  Vz' 
v„ 

x,  y,  z 
x',y  \z' 

Xj,  Xj 

Xi' 

xw 

ZS 

“5 

0£W 

A 

X 

V 

Ox,  Oy,  °z 

Subscripts 

ME 

MO 


Projections  of  V  onto  x',  y'  and  z'  axis  respectively,  ft/ sec 

x,  y,  and  z  components  of  V,  ft/ sec 

x,  y,  z  components  of  perturbation  velocity,  ft/sec 

f 

Free-stream  velocity,  ft/sec. 

Tunnel  coordinates,  rectangular  Cartesian,  in. 

LV  measurement  coordinates,  oblique  Cartesian,  in. 

Rectangular  Cartesian  coordinates  (i,  j  =  1,  2,  3),  in. 

Oblique  Cartesian  coordinates  (i=  1,  2,  3),  in. 

Chordwise  distance  from  wing  leading  edge,  in.  (Fig.  2c) 

Cartesian  coordinate  normal  to  local  surface,  in. 

Store  angle  of  attack,  deg 

Angle  of  incidence  between  wing  and  store,  deg 

Common  denominator  in  Eq.  (B-2) 

Wavelength 
Optical  frequency 

Standard  deviation  of  the  Vx,  Vy,  Vz  data  samples,  respectively,  ft/sec 

Mean 

Mode 


Abbreviations 


AP  Aperature 

BD  Beam  dump 


CF 


Color  filter 


DEC 


Dual  Bragg  cell 
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L 

LV 

M 

PMT 

R 

SBC 

W 


Single  Bragg  cell 
Window 
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